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Abstract 

We show that cosmic rays in external galaxies, groups and clusters rich in 
gas, with an average flux similar to that observed in the Milky Way, could 
have produced the observed extragalactic diffuse gamma radiation. 
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In addition to the galactic diffuse gamma radiation, which varies strongly with direction 
and can be explained by cosmic ray interactions in the galactic interstellar medium QXJ] , there 
appears to be an unaccounted diffuse component which is isotropic at least on a coarse scale 
and fits well at low photon energies to the extragalactic hard X-ray background radiation 
. These two features suggest an extragalactic origin of the isotropic diffuse gamma ray 
background radiation (GBR). Various unresolved extragalactic discrete and diffuse sources 
of gamma rays have been suggested |§ to explain its origin, its observed flux level ((5.5 ± 
1.3) xKT 5 and (1.3±1.3) x 1(T 5 7 cm" 2 s' 1 ster~ l above 35 MeV and 100 MeV, respectively) 
and its observed [f| spectral features (a power-law spectrum with a power index a = 2.4ta3 
between a few tens keV and a few hundreds Mev with a possible bump around a few MeV). 
The discrete sources included active galactic nuclei (AGN), Seyfert galaxies, BL Lac objects 
and normal galaxies while the diffuse sources included inverse Compton scattering of the 
microwave background photons in intergalactic space from cosmic ray electrons and more 
speculative sources such as matter-antimatter annihilation, cosmic strings, annihilation of 
super symmetric particles and decay of particles relics from the Big-Bang. However, all 
these explanations have been later questioned by observations. For instance, one of the 
more promising explanations of the GBR has been that it is the sum of gamma ray emission 
from unresolved AGN ||. When this explanation was suggested, only the relatively nearby 
quasar 3C 273 had been seen in high energy gamma rays M (E > 100 MeV). Since the 
launch of the Compton Gamma Ray Observatory (CGRO) in 1991, its Energetic Gamma 
Ray Experiment Telescope (EGRET) has detected 33 bright AGN in high energy gamma 
rays [0], all of which seem to belong to the "blazar" class. Many more AGN including 
blazars, which are both bright and relatively close and which were within the EGRET field 
of view, have not been detected in high energy gamma rays, indicating that not all of these 
objects are so luminous in gamma rays, or their emission is highly beamed or has a low duty 
cycle. The beaming hypothesis is further supported by other features of the sources such as 
superluminal velocity, beamed radio emission, and a power-law gamma ray spectrum (with 
a power index between 1.4 and 3.0, with values between 1.8 and 2 being most common). It 
is therefore difficult to see how the observed isotropic GBR can be produced by the highly 
beamed emission from AGN, many of which have a much harder spectrum than that of the 
GBR §. 

Most of the high energy gamma ray emission of our Milky Way (MW) galaxy can be 
explained by cosmic ray interactions with its interstellar gas jl]. While the total mass of 
gas in the MW is only about 10% of its stellar mass, i.e., M gas /M* « 0.10, recent X- 
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ray observations of clusters and groups indicate that they contain much larger fractions of 
gas ||. For instance, analyses of recent observations with the ROSAT X-ray telescope of 
the compact group HCG62 and the Coma cluster yielded M gas /M* rj 1.1 h~ 3 / 2 within a 
distance of 0.24/i" 1 Mpc from the center of HCG62 fTTJ, and M gas /M* « (9.0 ± 2.6)/i~ 3 / 2 
within a distance of 2.5/i _1 Mpc from the center of the Coma cluster |TT[, where H Q = 
100/i km s^Mpc" 1 is the Hubble constant. Detailed studies of the gas content of galaxies, 
groups and clusters indicate that for the whole Universe M gas ~ M*, i.e., M sos /M* = 
(3(M gas / M*)mw with (3 ~ 10. At higher redshifts a larger fraction of the baryons in the 
Universe are expected to reside in gas because cosmic evolution transforms gas into galaxies 
and stars. The gas content of the Universe at high redshift has been determined from the 
statistics of Lyman a absorption systems in quasar spectra measured with the large ground 
based telescopes and the Hubble Space Telescope. It is consistent with the present amount 
of visible matter in galaxies and clusters. In this letter we show that if the flux of high 
energy cosmic ray nuclei is either universal and equal to that observed in the MW, or if 
the MW cosmic ray flux represents well the average cosmic ray flux in clusters and groups 
of galaxies, then the observed extragalactic GBR could have been produced by cosmic ray 
interactions with gas in external galaxies, groups and clusters [|12[ . 

We first present two simple estimates of the extragalactic GBR produced by cosmic ray 
interactions in extragalactic gas which ignore cosmic evolution. Cosmic evolution effects 
(e.g., Hubble expansion, energy redshift, gas consumption by stellar formation, evolution 
of the cosmic ray flux) generally tend to increase the predicted GBR and will be included 
later. 

The main mechanisms by which cosmic rays produce high energy gamma - rays in the 
interstellar medium are bremsstrhalung and inverse Compton scattering from electrons and 
radiative decay of 7r°'s produced in inelastic collisions of cosmic ray nuclei (mainly protons) 
with the interstellar gas. Bremsstrahlung dominates gamma ray production at low energies 
while 7T° production dominates production of high energy gamma rays. Thus, for a universal 
cosmic ray flux, the emission is proportional to the mass of the gas (assuming a universal 
composition). The observed [Q total luminosity of the galaxy in gamma rays with energy 
above 100 MeV is L 7 = (2 ± 1) x 10 39 erg s~ 1 . The optical luminosity of the galaxy was 
estimated @ to be (2.3 ± 0.6) x 1O 1O L , while the average luminosity density of the Uni- 
verse was estimated [14]] to be (1.83 ± 0.35) x 10 8 h 2 L Q Mpc~ 3 . Consequently, cosmic ray 
interactions in extragalactic gas produce a GBR with an energy spectrum similar to that of 
the galactic diffuse gamma radiation and a flux level of, 
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It yields, e.g., F y (> 100 MeV^) ~ (5.0 ± 3.0) x 10~ 9 erg cm~ 2 s _1 ster' 1 , where we used 
the value, h ~ 0.80 ± 0.16, that was measured recently with the Hubble Space Telescope 



This flux level agrees, within the experimental uncertainties, with the flux level of the 
observed GBR. 

The large formal uncertainty in this estimate is due to the uncertain values of the total 
gamma ray and optical luminosities of the Milky Way and the value of the Hubble constant. 
If most of the baryons in the present Universe are still in the gas state then these uncertainties 
can be avoided by using an alternative estimate: 

F 7 (>E)^cH» 1 n b j H (>E). (2) 

For j H (E > 100 MeV) = (2.6 ± 0.24) x 10~ 26 7 cm~ 2 s' 1 ster' 1 (see ref. 0), where n b 
is the mean baryon density in the Universe and j H (> E) is the total emissivity in high 
energy gamma rays per baryon due to cosmic ray interactions. j H {> E) is known quite 
accurately, both theoretically and experimentally fl|], from cosmic ray interactions in the 
solar neighbourhood, j H (E > 100 MeV) = (1.8 ± 0.2) x 10 -26 7 cm~ 2 s~ x ster' 1 , while n b 
is known quite accurately from Big Bang nucleosynthesis (SBBN): The ratio of baryons 
to photons in the Universe, r\ = n b /n^ = (1.6 ± 0.1) x 10~ 10 , was deduced [|Tj| from the 
observed abundances of the light elements, H, D, 4 He and 7 Li extrapolated to zero age using 
Standard Big-Bang Nucleosynthesis (SBBN) theory. The microwave background radiation 
measured with the Cosmic Background Explorer (COBE) satellite is well fitted with a black 
body radiation of temperature [17]] T = 2.726 ± 0.017 K, yielding n 7 = 411 ± 8 cm -3 . 
Consequently, SBBN yields n b = (6.6 ± 0.4) x 10 _8 cm -3 and from Eq.2 it follows that 

F 7 (> 100 MeV) w (1.4 ± 0.3) x 10~ 5 7 cmT 2 s~ x ster' 1 , (3) 

which agrees with the flux level of the observed GBR. 

Note that we used a mean baryon density which is implied [|1^] by the recent mea- 
surements of the abundance of Deuterium in a high redshift low metallicity cloud [18| 
((1.9 — 2.5) x 10~ 4 by number) and also by the current best determinations [|l^] of the 



primordial abundance of 4 He (22.9 ± 0.5% by mass). However, if the primordial abundance 
of Deuterium is that observed in the local interstellar medium |2(J ((l-65l jg) x 10 -5 ) and if 
the current best determinations of the primordial abundances of 4 He and 7 Li underestimate 
their true values due to large systematic errors, then n b ~ 2.5 x 10~ 7 cm~ 3 and most of 
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the baryonic matter in the Universe is optically dark. Some of it, however, may reside in 
low density (X-ray dim) diffuse ionized gas. In that case Eq. 3 should be regarded as an 
estimated lower bound on the GBR. 

We will now consider cosmic evolution effects. In a Friedman-Robertson- Walker universe 
with a zero cosmological constant and a density Q in critical density units (i.e., Q = p/p c 
where p c = 3Hq/8tiG with G being Newton's constant of gravity), a population of gamma 
ray sources s with differential luminosities L S (E) oc E~ a , a number density at the present 
epoch n s and evolution functions f s (z) (which describe how the luminosities of sources s 
change with redshift z), generates a diffuse extragalactic flux at the present cosmic epoch, 
z — 0, which is given by 



In Eqs.1-4 we have neglected gamma ray absorption in the intergalactic space which becomes 
important only for E y > 500 GeV due to e + + e~ pair production in collision with photons of 
the intergalactic IR background radiation. Besides the gas content of the Universe the major 
unknown cosmological factor which affects Eq.4 is the evolution of the cosmic ray flux. We 
have investigated various plausible cosmic evolutions and found out that they can change 
the predicted flux level of the GBR as given by Eq.2 by less than a factor two. For instance, 
if high energy cosmic rays are universal and were injected into the intergalactic space during 
a short time period around redshift Zi (1 + Zi 3> 1), then their density at redshift z < zi was 
(1 + z) 3 times larger. Consequently, the emissivity of atoms was (1 + z) 3 times larger and 
the cosmic integration increases the GBR by a factor K pa (1 + Q/2)/(a + Q/2 — 2). The 
cosmic ray differential flux below 1000 GeV has a power index of a ~ 2.67, yielding K= 1.28 
for Q = 1 and K= 1.49 for Q = 0. 

Fig.l presents a comparison between our prediction of the cosmic ray produced extra- 
galactic GBR for the choice h = 0.8 and £1 = 0.15, and the observed extragalactic GBR. As 
explained, we estimate that the uncertainty in the absolute normalization of the flux level 
of the GBR due to uncertain cosmic evolution is less than a factor two. As demonstrated in 
Fig.l, an average cosmic ray flux in external, gas-rich, galaxies, groups and clusters, similar 
to that observed in the Milky Way can produce the observed extragalactic GBR. 

Our explanation of the GBR can be further tested in future observations, at least in two 
ways: 

a) At GeV energies gamma ray production is dominated by 7r° production. Because of 
Feynman scaling of pion production at GeV energies, the GBR should exhibit the same 
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power-law spectrum exhibited by the MW cosmic ray flux, namely, a ~ 2.67. Above 500 
GeV the power index should increase with increasing energy due to pair production off the 
extragalactic IR background photons. 

b) Nearby gas-rich clusters which were observed with the ROSAT X-ray telescope, such as 
the Coma cluster |TT[ (D m 70h~ 1 Mpc) and the Perseus cluster |2]J (D 55h~ 1 Mpc) may 
be detected by EGRET in high energy gamma rays. The expected gamma ray flux from the 
interaction of a MW like cosmic ray flux in a gas rich cluster at a distance D from Earth is 
given by, 

For the Coma cluster, M gas = (9.0±2.6) x lO 13 /r 5 / 2 M within a radius of 100 arcmin from its 
center |fLlf| . Consequently, its expected gamma ray flux level at Earth is, F 7 (> 100 MeV) ~ 
0.5 x 10~ 7 7 cm' 2 s _1 . For the Perseus cluster, M gas = (1.2 ± 0.3) x 1O 14 /i" 5/2 M within 
a radius of 100 arcmin from its center ||22|1 . Consequently, its expected gamma ray flux 
level at Earth is, F 7 (> 100 MeV) w 1.0 x 10~ 7 /i" 1/2 7 cm s . Such fluxes are marginally 
detectable by EGRET. However, a positive detection by EGRET of time independent gamma 
ray fluxes of these levels from the directions of the Coma and Perseus clusters will provide 
supportive evidence for the cosmic ray origin of the GBR. Conversely, a failure to detect 
gamma rays from the Coma and Perseus clusters by EGRET can be used to set a limit on 
the cosmic ray flux in these clusters, on the validity of cosmic ray universality and on the 
validity of the cosmic ray origin of the extragalactic GBR. 
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FIGURES 

FIG. 1. Comparison between the predicted (dotted line) extragalactic GBR produced by 
a universal MW-like cosmic ray flux in extragalactic gas (Eq.2) and the observed high 
energy GBR (dashed strip). 



9 
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